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The flow of motor proteins on a filamental track is modelled within the the framework of lattice 
driven diffusive systems. Motors, considered as hopping particles, perform a highly biased asym- 
metric exclusion process when bound to the filament. With a certain rate, they detach from the 
filament and execute unbiased random walk in the bulk which is considered as a closed cubic com- 
partment. Motors are injected (extracted) from the leftmost (rightmost) site of the filament located 
along the symmetry axis of the compartment. We explore the transport properties of this system 
and investigate the bulk-boundary interplay on the system stationary states. It is shown that the 
detachment rate notably affects the system properties. In particular and in contrast to ASEP, it is 
shown that the density profile of bound particles exhibit different types of non monotonic behaviours 
when the detachment rate varies. It is shown that in certain situations, the density profile of the 
filament consists of coexisting high and low regions. 
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I. INTRODUCTION 



Many transport problems in physics and related inter- 
disciplinary areas have been successfully modelled within 
the framework of lattice driven diffusive systems. The ki- 
netics of biopolymerization 1-3 , interface growth, 2 ' 4 and 
vehicular traffic flow 5-7 are among various examples of 
diverse applications of cellular automaton approach in 
describing out of equilibrium systems. Driven lattice 
gases with open boundaries are stochastic lattice mod- 
els of particles which hop preferentially in one direction 
and which are connected at their boundaries to parti- 
cle reservoirs of constant densities. As a result of this 
coupling, a stationary flow of particles is established in 
the system. A prototype example of non-equilibrium low 
dimensional systems is the asymmetric simple exclusion 
process (ASEP) where particles interact through hard- 
core repulsion 8 . ASEP has been applied to describe a 
variety of phenomena in physics and beyond. In one 
dimension, many interesting feature such as boundary- 
induced phase transitions, spontaneous symmetry break- 
ing and shock formation have been shown to exist 1 . 
These phenomena are associated to the boundary effects 
on the chain bulk and are generally termed boundary in- 
duced. Besides the boundaries, the bulk properties of 
the ID chain can be affected via coupling to a 3D bulk 
reservoir. A well known example is the Langmuir ki- 
netics (LK) 14 where particles can adsorb (desorb) if a 
bulk site is empty (occupied). Recently the interplay of 
both effects i.e., the boundary and the bulk reservoirs 
have received attention 15-18 . It has been verified that in 
certain situations where the adsorption/desopbtion rate 
properly scales with the system size, the density profile 
exhibits a non-monotonic structure in the bulk which is 



characterized by the coexistence of low and high den- 
sity regions separated by a shock front if the system size 
is high enough. In contrast to ASEP where shocks are 
not localized, in the presence of adsorption/desorption of 
particles, the system shows a more complex and richer 
behaviour. Specially, the shocks are localized. This 
behaviour can be verified by a continuous mean-field 
approach 15-17 . Along this line of study, we investi- 
gate the bulk-boundary interplay in a one dimensional 
chain. In our study, the homogeneous bulk reservoir is 
replaced by a more realistic structure. More specifically, 
the bound particles can detach from the one dimensional 
chain and execute unbiased random walk. The unbound 
particles are confined to a compartment with reflecting 
boundaries. It is one of our main objectives to see how 
the cooperative behaviours of the diffusive unbound par- 
ticles affects the results obtained from the static homo- 
geneous 3D bulk. We are interested to know if the shock 
localisation persists in the case where the bulk particles 
have dynamics. Our study is motivated from a biological 
application of ASEP. A considerable portion of intra-cell 
transport of biological cargos in eucaryotic cells is per- 
formed via the so-called motor proteins namely kincsin, 
myosin and dynein which move on filamentary tracks 
called microtubules and actins 19-26 . These bio motors 
perform a biased processive random walk on the filament 
and after some while detach from it due to thermal noise. 
Recently some physicists have attempted to formulated 
the problem described above within the lattice driven gas 
framework 27-31 . Besides investigating the diffusion prop- 
erties of a single molecular motor on a filament-like track, 
these investigations have been generalized to account for 
the mutual interaction of motors via hard core repulsion. 
The interaction was considered at its simplest form i.e., 
exclusion. It has been shown that anisotropy induced 
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by the filament can affect the diffusion properties of the 
bulk particles which arc unbound to the filament. On the 
other hand, the bulk diffusive motors can substantially 
influence the structure of the density and current pro- 
files on the filament. It has been recently been shown 32 
that in a compartment containing cytoskcletal filaments, 
hard core interaction between unbound motors, which 
are confined to move inside the compartment, and the 
filament gives rise to a localized shock structure in the 
density profile of the filament. The above geometrical 
arrangement i.e., a filament placed along the symmetry 
axis of a compartment are accessible to in vitro experi- 
ments and mimics the structure observed in cytosklcton. 
In particular, such systems serve as simple models for the 
description of motor-based transport in an axon 34 ' 35 . In 
this paper our main focus is on the role of detachment 
rate on the filament density profile. The investigation of 
the consequences of changing the detachment rate can be 
important not only from biological but also for theoreti- 
cal purposes. 



II. DESCRIPTION OF THE PROBLEM AND ITS 
FORMULATION 

Before proceeding further into the problem, let us 
define the model in more details. The motor proteins 
motion is assumed to take place in a compartment. For 
simplicity, we take the compartment shape a three di- 
mensional cube with dimensions Lx, Ly and Lz in the 
directions x, y and z respectively. We assume the fila- 
ment is located along the x axis of the tube. Further- 
more, the filament coordinates in y and z directions are 
taken as y — and z = k^-. We assume the particles 
(motors) enter onto the filament from the left boundary 
i.e., at the point (x, y, z) — (0, k^-). The particles can 
leave the compartment only via reaching the last point 
of the filament which is (x,y,z) — (Lx,^-,^-). The 
particles movement are governed according the following 
discrete-time stochastic dynamical rules. At each time 
step t, the state of the system is specified by the oc- 
cupation numbers of each site inside the compartment. 
If the site is full, its occupation number is one. If it is 
empty, the occupation number is zero. In addition, the 
motors are assumed to be point particles. Note that the 
integer coordinates i,j and k of the sites are restricted as 
< i < Lx; < j < Ly; < k < Lz. Now we explain 
the updating rules applied to all the particles inside the 
compartment. 

1) particle injection: 

In the beginning of the time step, we update the entry 
site of the filament: if the first site of the filament is 
empty, a particle will be injected to it with the proba- 
bility a. Furthermore, The number of entrants will be 
regarded as the system inflow. 



2) bound particle movement: 

We divide the particles inside the compartment into 
two groups: bound and unbound. If the particle is on 
the filament, we consider it as a bound one. If it is 
outside the filament, it is regard as an unbound particle. 
As discussed, the filament role is to provide a driven 
track on which the particles diffusive motion is highly 
directed. A particle on the filament can randomly move 
one step forward with probability p, one step backward 
with probability S, remains immobile at its position with 
probability 7 or leave the filament (detachment) to one of 
its non-filament neighbours with probability | 29 ' 30 . All 
these random movements will be successful if the target 
site is empty. In case the target site is already occupied 
by another motor, the attempted hop will be rejected. 
The probabilities should sum up to one: p+5+j+^- — 1. 
In three dimensions, there are four neighbouring non- 
filament sites adjacent to each filament site (i,-^-,^-) 
therefore the detachment probability equals ^ . 

3) unbound particle movement: 

An unbound particle is the one which is not located 
on the filament. These particles diffuse freely in the 
compartment space until they find a chance to attach to 
the filament. Once they attach to the filament, they be- 
come bound. The diffusion probabilities to the adjacent 
empty sites are equally taken to be i for all six spa- 
tial directions. For unbound particles, the compartment 
boundaries are assumed to be totally reflective. This can 
be implemented by putting immobile particles on the all 
boundary sites expect the entrance site (0, kj-, 4p) and 
the exit site (Lx, kj-, -^). 

4) particle extraction: 

Once a bound particle reaches the last filament site, 
it will be removed from the system (compartment) with 
the exit probability (3. The number of removed particles 
will be regarded as the system outflow. We recall that 
due to binding/unbinding of particles from the filament, 
in the steady state the current depends on the filament 
site number. Throughout this paper by current we mean 
the current of the last filament site. 

updating scheme: 

Before proceeding further, it would be illustrative to 
discuss, in some details, the type of our updating. Basi- 
cally, in the CA (cellular automata) models, the updating 
schemes are divided into two categories: particle- oriented 
and site- oriented. The particle oriented scheme can itself 
be divided into various methods. The most important 
one is parallel dynamics in which the updating rules are 
synchronously applied to all particles. The prototype ex- 
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ample is the vehicular traffic flow 5 . The next variant is 
the ordered sequential updating. In this scheme, one up- 
dates the particles in an ordered sequential manner say 
from left to right in one dimensional open systems. Fi- 
nally the third type is random sequential update. In this 
scheme, an update step is composed of many sub-update 
steps. In each sub-update step, we randomly chose one of 
the particles and update its position in accordance to the 
movement rules. In contrast to particle-oriented schemes, 
one has site-oriented ones. Site-oriented schemes, are 
generically divided into ordered and random sequential 
types. In the ordered scheme, one sweeps the lattice in a 
certain direction (prescribed by the method) and updates 
the position of those sites which are occupied by particles. 
In random sequential type, one randomly selects a site. If 
this site is occupied, then updating is done otherwise one 
continues by selecting other sites until an occupied site 
is found. The type of updating will affect some aspects 
of the problem but normally the phase structure of the 
model will not substantially changed under changing the 
updating scheme. In this paper we use particle-oriented 
ordered sequential scheme. More specifically, we label 
each particle upon entering into the compartment. At 
each time step, we update the particles' positions in the 
order of their labels. This resembles the site-ordered ran- 
dom sequential update. Our motivation to choose this 
update stems from our interest to evaluate the passage 
time of particles from entrance to removal. In the follow- 
ing sections we present our simulation results and discuss 
the phase structure and the other statistical properties of 
this system. 



III. SIMULATION RESULTS 

If we allow time to evolve for a long time we see the sys- 
tems reaches a non-equilibrium current-carrying steady 
state in which the statistical properties of the system 
does not show significant changes with time. To see this 
more explicitly, let us first define the quantities we are 
interested to compute: 

1) number of bound particles at the end of timestep t 
denoted by Nb(t). The subscript b denotes bound. Divid- 
ing Nb(t) by the filament size Lx gives the instantaneous 
global filament density at time t which we denote it by 
Pb(t). 

2) number of unbound particles at the end of timestep 
t denoted by N u b(t). Dividing N u b(t) by the compart- 
ment volume Lxx Lyx Lz~ Lx gives the unbound global 
particle density at time t which we denote it by p u b{t). 

3) filament current Jb(t) which is defined to be the 
number of particle exited from the last filament site by 
the end of t-th step divided by t. 



4) passage time of each particle: since we label the 
particles upon entering into the system, we can follow 
their trajectories from the time they enter the compart- 
ment until the time they leave it. Therefore we can 
specify each particle's passage time. p(t) is defined to be 
the average passage time of those particles which have 
already exited the compartment by the end of time t. 

5) centre of mass of the unbound particles denoted by 

< X u b > (t). This quantity is obtained by averaging 
over the ^-coordinate of all unbound particles at time t. 

6) centre of mass of the bound particles denoted by 

< Xb > (t) . This quantity is obtained by averaging over 
the ^-coordinate of all bound particles at time t. 

7) total number of particles in the system i.e., 
N u b{t) + Nb(t) which is named load. 

We note that in the steady state, the radial current 
of particles onto the filament should be equal to the cur- 
rent of particles leaving the filament. This is called radial 
equilibrium discussed in 29 ' 30 . In a mean-field approxima- 
tion, the radial equilibrium implies the following relation- 
ship between pb and p u b- 

|pt(l - Pub) = ^Pub(l - Pb) (1) 

6 6 

In the following we exhibit the stationary characteris- 
tics of the system. The system contains Lx = 100 sites. 
The perpendicular dimensions are Ly — Lz = 25 and 
5 = 7 = unless otherwise specified. In the first figure, 
we sketch the dependence of steady-state value of load in 
terms of a for three values of detachment rate e and exit 
rate /?. 





15000 
14000 
13000 
12000 
11000 




. 8 = 0.3 

< B = 0.3 

> 8 = 0.3 

1 8 = 0.7 

■ 8 = 0.7 

r p = 0.7 


E = 0.005 
8 = 0.010 
E = 0.015 
E = 0.005 

E = 0.010 A 
E = 0.015 




iady load 


10000 
9000 






J* / 




8000 
7000 












6000 
5000 
4000 
3000 
2000 
1000 






• '* * 







u o * 


0.2 


0.4 0.6 


0.8 1 










a 





Fig. 1: steady state value of load vs a. e and (3 are specified 
in the figure. 

We observe a phase-transition like behaviour near a c ~ 
0.3 and 0.7, for j3 — 0.3 and 0.7 respectively, where the 
steady load shows up a sharp increase. Moreover, in- 
creasing e gives rise to increment of the load. This is 
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expected since increasing c raises the detachment proba- 
bility hence the number of unbound particles is increased. 
To verify this transition, let us consider the dependence 
of the bound current Jj on a. 



see the phase transition is also manifested in the passage 
time. We note that transition is smooth for (3 = 0.7 and 
it resembles to a cross over. Let us next consider the 
density of bound motors p b as a function of a. 




Fig. 2: steady state current in terms of a for a system of 
size Lx = 100. The value of j3 and e are specified in the figure. 

for (3 = 0.3 the current linearly increases with respect 
to a and then saturates at an e-dependent value. The 
saturation value lies near 0.35. The smaller values of 
e correspond to higher currents. This can be explained 
on account of the fact that for small e motors prefer to 
remain bound to the filament track and continue their 
biased movement toward the exit point. Furthermore, 
we observe that the effect of increasing e is smoothing 
the transition. For larger value of (3 — 0.7, the overall 
picture remains the same. However the transition point 
and saturation value changes correspondingly. These be- 
haviour are analogous to ASEP. The next quantity is the 
passage time. We exhibit the steady state passage time 
vs a for some values of (3 and e. 



Fig. 4: steady state p b in terms of a for a system of size 
Lx = 100. The values of /3 and e are specified in the figure. 

Again, one observes phase transition behaviour from a 
low to a high density phase for (3 — 0.3 but a smooth 
crossover-like behaviour for (3 — 0.7. The results for 
(3 = 0.3 i.e., a first order low-to-high density transition is 
similar to that of ASEP. However when (3 — 0.7 one ob- 
serves slight differences to the behaviour of ASEP. Specif- 
ically we see deviations from linear dependence of p b vs a 
before reaching the saturation region. We also obtained 
the above diagram for Lx — 200 and 300. The results 
are very close to the case Lx = 100. We next explore 
the properties of the average position of the bound mo- 
tors < x > b . The behaviour of < x > b in terms of a is 
exhibited in the next graph. 




Fig. 3: steady state passage time in terms of a. Lx = 100 
and the values of (3 and e are specified in the figure. 



Fig. 5: steady state < x > b in terms of a for a system of 
size Lx — 100. The value of (3 and e are specified in the figure. 



The above graph is similar to the load graph. One can 



It is interesting to note that < x > b reaches its peak 
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at the transition point but decreases afterwards. Before 
transition, the system's load as well as pb are small and 
one does not expect jam formation hence < x >& should 
be around half of the filament length (here 50). Near 
transition point, jams are formed near the exit point and 
correspondingly < x >& increases over kf- = 50. The 
interesting point is that < x >b starts decreasing with 
further increase of a. This could be related to the bulk 
effects due to unbound particles. For (3 — 0.7 this effect 
is suppressed and the sharp peak is smeared out. The 
appearance of a peak in the centre of mass of bound par- 
ticles signifies an important characteristics. The peak 
tells us that around the value a ~ 0.35 an inhomoge- 
neous density structure forms throughout the filament. 
This point will become more clear by examining the fila- 
ment density profile. To get a deeper insight, it would be 
useful to consider density profile on the filament. Figure 
(6) shows such profiles. 




site number 



Fig. 6: density profile on the filament for various values of 
a for a system of size Lx = 100. j3 = 0.3 and e = 0.005 . 



For small values of e, we observe analogous behaviour 
to ASEP. One has an exponential increase near the right 
boundary in the low density phase. In marked contrast 
to ASEP, we see that for a equal to 0.35, there is a coex- 
istence of low and high density regions. In other words, 
there is a localized shock structure. This has been ear- 
lier reported when the bulk was modelled by constant 
attachment/detachment rates. 15-17 . Due to finite size of 
the system, the domain wall is not so sharp. Upon in- 
creasing a, the density again becomes nearly constant in 
the bulk with an exponential density decrease near the 
left boundary. To see the effect of varying e, in the fol- 
lowing figure, we have depicted the density profile for the 
same values of a but with e = 0.015. 
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Fig. 7: density profile on the filament for various values of 
a for a system of size Lx = 100, (3 = 0.3 and e = 0.015. 

Analogous to the case e = 0.005, we again have a lo- 
calized shock structure at the same a but this time the 
position of the domain wall is shifted towards left. Fur- 
thermore, one can recognize a linear profile, with a small 
slope, in the high density regime. Although the slope 
is small but it is yet a distinguishing aspect. This can 
be related to formation of spontaneous jam or moving 
shocks and is a new feature of the problem which is ab- 
sent in ASEP. In fact, the interaction of one dimensional 
filament with the unbound bulk particles leads to this 
weakly linear density profile. Now let us investigate the 
model properties when the exit rate (3 is changed. The 
following figure illustrates the behaviour of steady load 
vs [3 for a = 0.3 and 0.7 each for three values of e. 




Fig. 8: steady state load vs (3 for a system of size Lx = 100 
for various values of a and e. 

For a — 0.3 one observes a transition-like behaviour 
while for larger a = 0.7 the abrupt behaviour is changed 
into a smooth one. The effect of increasing e is again 
smoothing the transition behaviour. The passage time 
diagram is similar to the load diagram and we do not 
show it. More interesting is the steady current behaviour. 
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Fig. 9: steady state current in terms of for a system of 
size Lx = 100. The value of a and e are the same as in the 
previous figure. 

Analogous to ASEP, the current saturates at a in- 
dependent value. The effect of e is notable only for high 
a, corresponding to high density phase, where its increase 
gives rise to a current decline as expected. Besides this 
feature, for high a we see that varying e gives rise to 
smoothing of the changes in the current. Next we look 
at the global density on the filament. The following graph 
depicts this behaviour. 
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Fig. 10: steady state pb in terms of for a system of size 
Lx = 100. The value of a and e are specified in the figure. 

Similar to ASEP, one observes a first order high-to-low 
density transition for a = 0.3 and a smoothly decreasing 
behaviour for a = 0.7 until it reaches a constant value. 
The transition occurs at ~ 0.25. We also examined the 
case Lx — 200 and 300. The results do not show signifi- 
cant changes to the case Lx =100. This shows that the 
behaviour is not due to finite size effects but is related to 
the bulk influence. We have also evaluated the < x >& 
i.e., the average position of bound particles vs 0. Figure 
(11) shows this dependence. 



Fig. 11: steady state < x > b in terms of for a system of 
size Lx = 100. The value of a and e are specified in the figure. 



Similar to figure (5), we see a sharp increase of < x >b 
to a maximum value ~ 0.6Lx at the transition point. It 
decreases to half of the filament for larger values of 0. 
The position of maximum crucially depends on a and 
e. For larger values of a, the maximum point shifts to- 
wards larger 0. Moreover, the maximum value of < x >b 
reduces. The reason is that for larger a, the system 
contains more particles and can maintain high current. 
This leads to a more homogenised density profile. Lastly, 
we depict the density profiles for some values of with 
e = 0.01 and 0.015 in the following figures. In comparison 
to ASEP, we see entirely different profiles. Generally, the 
profile is increased and then starts declining in a linear 
manner. 
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Figs. 12: density profile on the filament for various values 
of for a system of size Lx = 100 and a — 0.7. e = 0.01. 
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Figs. 13: density profile on the filament for various values 
of (3 for a system of size Lx — 100 and a = 0.7. e = 0.015. 

The effect of increasing (3 is two fold. First, smoothing 
the changes between increasing and decreasing regions 
of the density profile and second, extending the region 
of increasing profile. In comparison to the top figure 
where e = 0.01, for a larger value of e = 0.015 we see 
that the size of region where the profile is increasing, de- 
creases. Increasing e is more effective when the exit rate 
(3 is large. For instance, when (3 = 0.63, the increasing 
part of the density profile decreases in comparison to the 
case e — 0.005. 

To see the behaviour of the density profile for larger 
systems size, we also obtained the profile for the same 
parameters as in the above figures but this time for 
Lx = 300. The following figure depicts the behaviour: 




Fig. 14: density profile on the filament for various values 
of j3 for a system of size Lx — 300, a = 0.7 and e = 0.015. 

One observes that for large (3 a linear profile forms in 
the bulk. The effect of the left boundary on the profile 
is the formation of an increasing density layer. The loca- 
tion of the maximum point moves away the left boundary 
for higher (3. One also observes a concave density profile 
structure for sufficiently small (3. This aspect was absent 



in Lx = 100 where not only the concavity did not exist 
but also for (3 — 0.77 there existed a concave profile. In 
order to have a deeper insight to the problem, we looked 
in more details into the effect of e on the density profile. 
In the following figures, we show the density profile for 
some values of e . a and (3 are fixed. We first consider 
the case a = (3 = 0.7. 
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Fig. 15: density profile on the filament for various values 
of e for a system of size Lx = 200, Ly = Lz — 25, a = 0.7 
and (3 = 0.7. 

Except, a short region where the profile is increasing, 
we again observe a linearly decreasing profile. The over- 
all effect of increasing e is to reduce the density value 
throughout the chain. The reason for such decrease in 
the local densities could possibly be explained on account 
that in the present case a = (3 the system can maintain a 
large current therefore increasing e would effectively re- 
duces the number of bound particles and therefore the 
local density on the filament decreases. The next figure 
exhibits the density profile for the case a = (3 = 0.3 . 
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Fig. 16: density profile on the filament for various values 
of e for a system of size Lx = 200, Ly = Lz — 25, a = 0.3 
and (3 = 0.3. 
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Here we observe a constant profile in the bulk accom- 
modated with sharp increasing layer at the right bound- 
ary. This is also in contrast to ASEP where once expects 
a linear profile along the a = (3 < \ line. We conclude 
that the the detachment rate e can give rise to signifi- 
cantly different features which are absent in ASEP. This 
conclusion have been earlier obtained for similar mod- 
els without taking into account the dynamic of unbound 
particles 15 " 18 ' 32 . Our final graph depicts the density pro- 
file for some values of larger e. a and (3 are fixed but 
unequal. 




Fig. 17: density profile for various values of e for a system 
of size Lx = 200, Ly = Lz = 25, a = 0.3 and (3 = 0.7. 

The above graph exhibits the profile in a low conges- 
tion state corresponding to a = 0.3 and [3 — 0.7. As 
observed, the effect of e is more notable for higher val- 
ues. For small e, the effect is an upward shift in the 
density profile. Larger e not only shifts the density pro- 
file towards a higher value but also creates a slope in it. 
Finally, to see the effect of lateral dimensions of the com- 
partment i.e., Ly and Lz, we performed extensive Monte 
Carlo simulations for Ly = Lz = 12 up to 70. Our re- 
sults showed no significant changes for the case where e 
is small. 



have studied the effect of varying the detachment rate e 
on the transport characteristics of the model. Although 
the model behaviours resemble that of ASEP in some 
ranges of parameters, especially when the unbinding rate 
is sufficiently low, Our investigation has illustrated that 
for considerable high values of the unbinding rate, one ob- 
serves remarkable distinct features with respect to ASEP. 
Specifically, the density profile shows linear and localized 
behaviour in some ranges of parameters. This implies 
that the density shock structure of the models should 
be entirely different than ASEP. This aspect can be of 
relevance to biological applications. Besides biological 
motivation, our study shed more lights onto the prob- 
lem of bulk-boundary interplay in the transport char- 
acteristics of low dimensional non-equilibrium systems 
which recently has attracted the attention of researchers 
in the field. Analogous to the cases where the detach- 
ment/attachment rates are constant 15-17 , we see the co- 
existence of high and low density regions which is termed 
shock localization. Moreover, we also observe other types 
of density structures such as linear density profiles with 
distinctive boundary behaviours. This shows that incor- 
porating dynamics for the bulk particles, can drastically 
affect the phase structure of the system. It would be 
illustrative to find out the phase diagram structure of 
the model. This is our next step. Finally we recall that 
our model does not incorporate the internal degrees of 
freedom for motors. In practice, these motors resemble 
a two-headed creatures. Taking this fact into a quanti- 
tative description 36 together with chemical state of the 
motors are two crucial steps which should be considered 
in further progresses 37 . Work along this line is in prepa- 
ration and will be reported elsewhere. 
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IV. SUMMARY AND CONCLUDING REMARKS 



Let us now summarize what has been done in this pa- 
per. We have developed a discrete time cellular automata 
model for the description of motor protein traffic through 
a filament-like track. Motors execute totally asymmetric 
random walk when bound onto the track. They ran- 
domly detach from the filament and perform unbiased 
random walk inside a closed compartment with reflect- 
ing boundary conditions. A similar problem has been 
earlier considered in 29 . Our boundary condition and up- 
dating schemes is different to that in 29 . Our main interest 
has been the interplay of the bulk and the boundary on 
the stationary properties of the system. In particular we 
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